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Abstract. We survey some aspects of the theory of the integrated 
density of states (IDS) of random Schrodinger operators. The first part 
motivates the problem and introduces the relevant models as well as 
quantities of interest. The proof of the existence of this interesting 
quantity, the IDS, is discussed in the second section. One central topic of 
this survey is the asymptotic behavior of the integrated density of states 
at the boundary of the spectrum. In particular, we are interested in 
Lifshitz tails and the occurrence of a classical and a quantum regime. In 
the last section we discuss regularity properties of the IDS. Our emphasis 
is on the discussion of fundamental problems and central ideas to handle 
them. Finally, we discuss further developments and problems of current 
research. 
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1. Introduction 

1.1. Physical Background: Models and Notation. The time evo- 
lution of a quantum mechanical state ip is obtained from the time dependent 
Schrodinger equation 

d 

i-^^P = Hrp. (1.1) 
Another important differential equation is the heat or diffusion equation 

^il>=-Hi/>. (1.2) 

In order to investigate either of these equations, it is extremely useful to 
know as much as possible about the operator H and its spectrum. In general 
the Schrodinger operator H is of the form 

H = H + V. 

The free operator Hq represents the kinetic energy of the particle. In the 
absence of magnetic fields, it is given by the Laplacian 

d FP- 
# n = _A = - V 



i dx u 2 

The potential V encoding the forces F(x) = —W(x) is acting as a multi- 
plication operator in the Hilbert space I? 1 



Occasionally, we will replace Hq with the operator Hq(B) which contains 
a homogeneous magnetic field. For d = 2 Hq(B), B > 0, is given by 



d 1 \ 2 / d 1 



In contrast to Hq(B = 0), the spectrum of Hq{B) for B > is a countable 
set a(H (B)) = {(2n + l)B;n E N}. The energies E n = (2n + 1) are 
eigenvalues of Hq{B) of infinite multiplicity, called the Landau levels. 

1.2. Random Potentials. First-order physical modeling often assumes 
an ideal background, e.g., a homogeneous material without any impurities. 
For example, in ideal crystals the atoms or nuclei are supposed to be dis- 
tributed on a periodic lattice (say the lattice Z rf ) in a completely regular 
way. We assume that a particle (electron) at the point x S M rf feels a poten- 
tial of the form q f(x — i) due to an atom (or ion or nucleus) located at the 
point i G Z d . We call the function / the single site potential. The coupling 
constant q represents the charge of the particle at the lattice point i. So, in 
a regular crystal, the particle is exposed to a total potential 



V(x) - qf(x- i). 



(1.3) 
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The potential V in (jOj) is periodic with respect to the lattice Z d , i.e., V(x — 
i) = V(x) for all x 6 M d and i G Z d . The mathematical theory of Schrodinger 
operators with periodic potentials is well developed (see, e.g., \35\ \112\ 
I15()j ). It is based on a thorough analysis of the symmetry properties of 
periodic operators. The spectrum of periodic Schrodinger operators has a 
band structure and the spectrum is purely absolutely continuous, i.e., 

oo 

a{H) = (J [ with a n < b n < a n+ i, (1.4) 

n=0 

(H) = 0, a(H) = a &c (H). (1.5) 

The real world is not ideal. Solids occur in nature in various forms. Some- 
times they are (almost) totally ordered, sometimes they are more or less 
completely disordered. For a mathematical modelling of disordered solids 
two ingredients are essential: the spatial homogeneity in the mean and the 
disappearance of long range correlations. In full generality these properties 
are studied in the theory of ergodic operators. This class of operators con- 
tains, e.g., quasicrystals and can also be related to random matrices. In this 
publication we will be concerned with random Schrodinger operators. 

The most popular and best understood model of a disordered solid is 
the alloy-type Schrodinger operator. It models a mixture of different atoms 
located at lattice positions. The type of atom at the lattice point i is assumed 
to be random. These particles are represented by randomly distributed 
coupling constants q± encoding the different charges. The total potential is 
given by 

V u (x) = Y, *(")/(*-«)■ (1-6) 

When talking about the alloy-type model, we will mean (jTJj) with the fol- 
lowing assumptions: 

(1) The single site potential / is bounded, non-negative and strictly 
positive on an open set. 

(2) / satisfies fix) < C (l + | x|)-( d+£ ) for some C and e > 0. 

(3) The random variables qi are independent and identically distributed 
random variables on a probability space (O, P). 

(4) The common probability distribution of the qi is denoted by -Po- 
rts support supp Pq is compact and contains at least two points. 

Assumption can be considerably relaxed. For example, for most of the 
following results one may allow local singularities for /. By assuming Q 
we avoid technical difficulties which may obscure the main argument. More 
details on weaker conditions can be found in the papers cited. Assumption 
((2j) ensures that the sum in (|1.6|) is convergent. The compactness of suppPo 
is convenient but not always necessary and in some especially marked sit- 
uations we consider also unbounded single site distributions. However, for 
many of our results we need that the qi (and hence suppPo) are bounded 
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below. The physical model suggests that suppPo consists of finitely many 
points only (the possible charges). However, for many mathematical results 
it is necessary (or at least convenient) to suppose that Pq has a (bounded) 
density g, i.e., Pq = g(X) d\. 

One might argue that such an assumption is acceptable as a purely tech- 
nical one. On the other hand one could argue the problem is not understood 
as long as it is impossible to handle the physically relevant case of finitely 
many values. 

A simplified version of the alloy-type Schrodinger operator above is the 
(discrete) Anderson model. Here the Hilbert space is the sequence space 
£ 2 (7j d ) instead of L 2 (R d ) and the free operator H$ is replaced by the discrete 
analogue of the Laplacian. This is the finite-difference operator 

(hou)(n) = — (u(m)—u(n)). (1-7) 

| TTi— n|=l 

Above we set |n| = Yli=i l n «l on The potential is a multiplication 
operator V = on i 2 (7L d ) with V^ii) independent, identically distributed, 
and the total Hamiltonian is defined by h = ho + V u . We will call this setting 
the "discrete case" in contrast to Schrodinger operators on L 2 (R d ) which we 
refer to as the "continuous case." 

The most frequently used and most important approach to model amor- 
phous material like glass or rubber defines the random locations of the atoms 
by a Poisson random measure. This random point measure can be character- 
ized by the number ua = ^(A) of random points in the set A. We assume 
that the random variables ha and ub are independent for disjoint (mea- 

surable) sets A, B and P (ua = k) = -A- e~ 1 (\A\ denotes the Lebesgue 
measure of A). With this notation we may write the Poisson potential for 
an amorphous solid as 

V u (x)= ( qfix-^d^irj). (1.8) 

To model thin disordered layers, we also consider random potentials which 
are concentrated along a hypersurface in M. d (or Z d ). For example, we are 
going to consider "surface" alloy potentials. To define such a potential let 
us write R d = R dl x R d2 then 

V u) (x 1 ,x 2 )= ^2 1ii( UJ )f(xi-h,%2) 

is a random potential which is concentrated along the hypersurface R dl in 
R d . In addition, there may be a random or periodic background potential 
on R d . 

Most of the theorems we are going to discuss can be proved for rather 
general single site potentials / and probability distributions Pq of the qi. 
For example, most of the time we can allow some local singularities of /. To 



6 



WERNER KIRSCH AND BERND METZGER 



simplify the following discussion, we will assume in this paper the conditions 
defined in the context of the alloy-type potential. 

The above random operators are examples of "ergodic operators." This 
class of operators includes not only most random operators but also periodic 
and almost periodic operators. Most of the results of Section |2] and part of 
Section H] can be shown for general ergodic operators. We refer to [148\ IT5^ 
12 7^ 11191 1166) and |68j for a discussion of this general context. 

1.3. The Concept of the Integrated Density of States. The (in- 
tegrated) density of states is a concept of fundamental importance in con- 
densed matter physics. It measures the "number of energy levels per unit 
volume" near (resp. below) a given energy. 

Typical systems arising in solid state physics have periodic or ergodic 
potentials. Consequently, the spectrum of the corresponding Hamiltonian 
is not discrete. Therefore, we cannot just count the eigenvalues below E 
or within an interval \E\,E2\. On the other hand, the number of electrons 
in such a system, which extends to infinity, ought to be infinite. For these 
two reasons, the Pauli exclusion principle does not make immediate sense. 
(How do we distribute infinitely many electrons on a continuum of spectral 
energies?) 

However, there may be a chance to make sense out of the Pauli principle 
by first restricting the system to a finite volume A. Inside A there should 
be only finitely many electrons, in fact, we may assume that the number of 
electrons in a given A is proportional to the volume |A| of this set. 

If P is a finite-dimensional orthogonal projection, then tr(P) is the 
dimension of its range. If P(~oo,E] i s the spectral projection of a ran- 
dom Schrddinger operator (which as a rule has infinite-dimensional range) 
and if A^ is a cube of side length L around the origin, then we may call 
tr(xA L P{-oo,E]) the restriction of P(-oo,E] to the cube A^. Above xa is 
the characteristic function of the set A. Thus, we may try to define the 
integrated density of states as 

N(E) = lim -L tr { X A L P(-oo,E]) • (1.9) 

Li— >oo -li-Zj 

Of course, we have to prove that the limit in (|1.9|) does exist and is not 
trivial. We will deal with these questions in Section 2. 

There is another way to define the integrated density of states which 
turns out to be equivalent to ()1.9|) : We restrict the operator H w to the 
Hilbert space L 2 (A). To obtain a self adjoint operator we have to impose 
boundary conditions at <9A, e.g., Dirichlet or Neumann boundary conditions. 
We call the corresponding operators and H^, respectively. These op- 
erators have compact resolvents, i.e., their spectra are purely discrete. We 
denote by 

Ei(fff) < E 2 (H°) < E 3 (H°) . . . (1.10) 
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the eigenvalues of (and analogously for Hj^) in increasing order, where 
eigenvalues are repeated according to their multiplicity. The eigenvalue 
counting function of an operator A with purely discrete spectrum is defined 
by 

N(A, E) = #{n\ E n {A) < E} = tr {P^e] (A)) . (1.11) 
Analogously to (|1.9|) . we can therefore define 

N D (E)= km J-.N{Hg,E) 



L^oo |A 



L 



= lim -L tr^.^tffj) (1.12) 

and similarly for Neumann boundary conditions, 

N N {E) = lim j^- NtH^E) 

= lim -L tr(P ( _^(fl£)). (1.13) 

This procedure to define the integrated density of states makes sense only 
if N, N D and N N all exist and agree. 

This is, indeed, the case. We will see in the sequel that each of these 
definitions has its own technical advantage. The integrated density of states 
N is basic for studying the physical (in particular the thermodynamical) 
properties of disordered systems. ^From a mathematical point of view, the 
properties of N are interesting in their own respect. Moreover, properties 
of N constitute an essential input to prove localization properties of the 
system. 

It is the aim of this review to discuss some of the problems and results 
connected with the integrated density of states. In Section 2 we sketch 
the proof of the existence of the integrated density of states and discuss 
some fundamental questions concerning the probabilistic and the functional 
analytic approach. In Section 3 we study the behavior of the integrated 
density of states at the boundary of the spectrum. In the last section we 
discuss some basic ideas concerning the regularity of the integrated density 
of states. 

Acknowledgements. It is a pleasure to thank many colleagues for fruitful 
collaborations and stimulating discussions on the subject. There are too 
many to name them here. We also would like to stress the fact that the 
selection of topics within our subject and the way of presenting them is due 
to our very personal preferences. We have most certainly left out important 
topics and works. This is to be blamed on our ignorance and the limitation 
of time and space. 

We would like to thank Jessica Langner and Riccardo Catalano for their 
skillful typing and careful proofreading of the manuscript. 
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2. The Density of States Measure: Existence 

2.1. Introduction and Historical Remarks. The first existence proofs 
for the integrated density of states go back at least to Pastur. See |145j for 
an early review of the subject. 

There are a couple of methods to prove the existence of the integrated 
density of states. One of them, invented and used by Pastur, is based on the 
Laplace transform of the integrated density of states and of its approximants. 
For this method, one proves the convergence of the Laplace transform and 
uses the fact that convergence of the Laplace transform implies the vague 
convergence of the underlying measures. 

To prove the convergence of the Laplace transforms, it is useful to express 
the Laplace transform of the finite-volume quantities using the Feynman- 
Kac-representation of the Schrodinger semigroup e~ tH . Feynman-Kac and 
Laplace transform methods are also used to prove the equivalence of the 
definitions of the integrated density of states 1)1.9)1 and ()1.12[) . 1)1.13(1 with 
either Neumann or Dirichlet (or more general) boundary conditions (see, 
e.g., Pastur [TI5] or [THESE!]) 

The definition of the integrated density of states via 1)1.9)1 was used by Avron 
and Simon in the context of almost periodic potentials j^j- They also proved 
that the spectrum of the operator coincides with the growth points of the 
integrated density of states. In Section [2.21 we will follow this approach to 
prove the existence of the integrated density of states. 

One of the virtues of the definition of the integrated density of states via 
boundary conditions 1)1.12)1 and 1)1.13)1 is the fact that they allow lower and 
upper bounds of N "for free." In fact, one way to examine the behavior of 
N at the bottom (or top = oo) of the spectrum is based on this approach. 
We will discuss this approach in Section 12.31 and the estimates based on it 
in Section 3. 

As a rule, quantitative estimates on the effect of introducing boundary 
conditions are hard to obtain. For example, if one investigates the behav- 
ior of N at internal spectral edges it seems extremely difficult to control 
the perturbation of eigenvalues due to boundary conditions. Klopp |97| 
proposed an approximation of the random potential by periodic ones with 
growing period. This way we lose monotonicity which is at the heart of the 
Neumann-Dirichlet approach. Instead one can prove that the approxima- 
tion is exponentially fast thus gaining good estimates of the remainder. 

Finally, we would like to mention that one can also define the integrated 
density of states via Krein's spectral shift function. This reasoning is well 
known in scattering theory (see, e.g., )9"|, 1180) ). In connection with random 
Schrodinger operators, the spectral shift function was first used by Simon 
)159j to investigate spectral averaging. Kostrykin and Schrader )106j ap- 
plied this technique to prove the existence of the integrated density of states 
and the density of surface states. This method turns out to be useful also 
to investigate regularity properties of the integrated density of states )107j . 
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The results of this section are true not only for the specific random 
potentials discussed in Section Q but rather for general ergodic operators. 
In fact, the proofs carry over to this general setting in most cases. We refer 
to the survey |68j for details. 

2.2. The Existence of the Integrated Density of States. In this 
section we prove the existence of the integrated density of states as defined 
in 1)1.9)1 . To do so, we need little more than Birkhoff's ergodic theorem 
(see, e.g., |113j ). Below, as in the rest of this paper, we denote by E the 
expectation with respect to the probability measure P. 

Proposition 2.1. If ip is a bounded measurable function of compact 
support, then 

lim -L tT(ip(H u )xA L ) =E(tr( X Ax¥>(#u,)XA 1 )) (2.1) 
for W-almost all lo. 

Proof. Define = tr (^(H^XAxii))- d is an ergodic sequence of ran- 
dom variables. Hence, Birkhoff's ergodic theorem implies that 

-i- tT(ip(H„) X A L ) = iTT £ & 
|Al| |A L | 

converges to its expectation value. □ 

The right-hand side of 1)2.1(1 as well as | | — 1 tr ( ^{Hj) XA L ) are posi- 
tive linear functionals on the bounded, continuous functions. They define 
positive measures v and vl by 

^ <p(\)du(\) = E( ( tr ( XAl tp(H u ) X A 1 )) ) 
and 

1 



ip(X) du L {\) = ——\x{(p{Hj)xh L ). 

Equation (|2.1|) suggests that the measures vl might converge to the limit 
measure u as L — > oo in the sense of vague convergence of measures. The 
problem is (|2.1|) holds only for fixed ip on a set £1^ of full probability; respec- 
tively 1)2.1)1 holds for all for w £ f^^O,^. However, this is an uncountable 
intersection of sets of probability one. The problem is solved by approxi- 
mating Co(R) by a countable, dense subset. 

Theorem 2.2. The measures vl converge vaguely to the measure v P- 
almost surely, i.e., there is a set of probability one, such that 

p{\)dv L {\) J <p{\)dv(\) (2.2) 

for all ip £ Cq(M), the set of continuous functions with compact support, and 
all u £ Qq. 
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Definition 2.3. The non-random probability measure v is called the 
density of states measure. The distribution function N of v, defined by 

N(E)=v(]-oo,E]), 

is known as the integrated density of states. 

Using Theorem 12.21 it is not hard to see: 

Proposition 2.4 ([3]). supp(z^) = S [= cr{H u ) a.s.]. 



2.3. Existence via Dirichlet Neumann-Bracketing. Our first ap- 
proach to define the density of states measure was based on the additivity 
of tr((^(i? a; )^A L ) an d the ergodic theorem by BirkhofF. This very naturally 
fits in the concept of self- averaged quantities from physics. 

However, for some part of the further analysis, an alternative approach — 
the Dirichlet-Neumann bracketing — is more suitable. Let (-Hoj)a and (H^)® 
be the restrictions of H u to L 2 (A) with Neumann and Dirichlet boundary 
conditions. See, e.g., [150] for an appropriate definition of these boundary 
conditions via quadratic forms. Furthermore, we define (for X = N or D, and E 6 
R) 

N*(E) := N((H U )%,E) = tr( X( ^ E] (H^)). (2.3) 
Our aim is to consider the limits 

N X (E) = lim -J-jV£(£). (2.4) 

The quantities and as defined in (|2.3|) are distribution functions of 
point measures vff and concentrated in the eigenvalues of and Hj^ , 
i.e., 

N£(E) = v£((-oo,E}). (2.5) 

The convergence in (f2.4|) is meant as the vague convergence of the corre- 
sponding measures or, what is the same, as the pointwise convergence of the 
distribution function at all continuity points of the limit. 

Let us first look at m (E) . The random field is not additive 
in A, so that we cannot use BirkhofF's ergodic theorem. However, is 
superadditive, in the sense that (E) > (E) + N® 2 (E) whenever 
A = Ai U A 2 with (Ai)° n (A 2 )° = 0. (M° denotes the interior of the set M.) 
Similarly, is subadditive, i.e., —N N is superadditive. 

Theorem 2.5. is superadditive and is subadditive. More pre- 
cisely, if A = Ai U A 2 and (Ai)° n (A 2 )° = then 

N^(E)+Ng(E) < N°(E) < <(£) < <(£) +Ng(E). 
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Fortunately, there are sub- and superadditive versions of the ergodic 
theorem, going back at least to Kingman |66j . The situation here is ideal 
for the superadditive ergodic theorem by Akcoglu and Krengel [2]. Indeed, 
one can prove that (for fixed E) the processes N A {E) and N A (E) are su- 
peradditive and subadditive random fields in the sense of [2] respectively 
(see |74l 1111) )- This yields the following result. 

Theorem 2.6 ([Sj). The limits 



and 



N D (E)=lka —N(H° A ,E) 

L->oo | Ai| 



N N (E)= lim —N(H» A ,E) 



exist F '-almost surely. Moreover, 

N D (E) = sup -i- E(N(H2 Al ,E)) 

L | ^ *-L 1 

N N (E) = M JL E(N(H% L ,E)). 

The functions N x are increasing functions. However, it is not clear 
whether they are right continuous. To obtain distribution functions, we 
define N x by making the N x right continuous 

N D (E) = inf N D (E') (2.6) 

E'>E 

N N (E) = inf N N (E'). (2.7) 

E>>E 

Note that N x and N x disagree at most on a countable set. Since N D are 
N N are distribution functions, they define measures by 

u D ((a, b}) = N D (b) - N D {a) (2.8) 

u N ((a,b}) =N N (b)-N N (a). (2.9) 

^From Theorem 12. 61 we obtain the following corollary, which we will use 
to investigate the asymptotic behavior of the integrated density of states 
(e.g., for small E). 

Corollary 2.7. For any A, 
^-E(N(H^ L ,E))< N D (E) < N N (E) < ±-E(N(H^ L ,E)). 

Our physical intuition would lead to the hope that N(E) = N D (E) = 
N N (E) since, after all, the introduction of boundary conditions was a math- 
ematical artifact that should not play any role for the final physically mean- 
ingful quantity. This is, in fact, true under fairly weak conditions (see 
|145l 1741 134j and references given there) . 
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Theorem 2.8. The distribution functions N(E), N (E) and N (E) 
agree. 

Theorem 12 . 81 follows from Theorem l2.1Ul in the next section. An alterna- 
tive proof for the Anderson model can be found in the review |7U| . Theorem 
|2~H1 implies a fortiori that the quantities -n^-r Nf? L (E) and -r±j Ng L (E) 
converge to the same limit, except for a countable set of energies E. The 
exceptional points, if any, are the discontinuity points of N. We will discuss 
continuity (and, more generally, regularity) properties of N in Section 4. 

2.4. A Feynman Kac Representation for N. In this section we 
will consider the Laplace transform of the integrated density of states (both 
N(E) and N X (E)). The Laplace transform of a measure v with distribution 
function F is defined by 

u(t) = F(t) := J e~ xt dv{\) = j e~ xt dF{\). (2.10) 

There is a very useful representation of the Laplace transform of N 
(and of N x ) via the Feynman-Kac formula. Using this representation one 
can show that N and N x are, indeed, the same quantities. Moreover, the 
Feynman-Kac-representation of N is very useful to compute the asymptotic 
behavior of N for small or large energies. 

The key ingredients of the representation formula for N(t) are the Brow- 
nian motion, the Brownian bridge and the Feynman-Kac- formula. For ma- 
terial about these concepts, we refer to Reed-Simon [1491 1150j and Simon 

HM|. 

By Pq^ we denote the measure underlying a Brownian bridge starting 

in the point x G M. d at time and ending at time t in the point y. Eq^. 

denotes integration over Pq^- A Brownian bridge is a Brownian motion 

b conditioned on b(t) = y. Note that Pg^. is not a probability measure. 

Fq V x has total mass p(t, x, y) where p denotes the probability kernel of the 
Brownian motion. 

Theorem 2.9 (Feynman-Kac formula). If V £ L[ oc unif (]R d ) for p = 2 

ifd<3,p> d/2 ifd>3, then e~ tH has a jointly continuous integral kernel 
given by 

e~ tH (x, y) = 4'£(e-/oM*M) ( 2 .H) 
= J e-^ v ^ ds dF^ x {b). 

The integration here is over paths b{.) G C([0,i]). 

We remind the reader that we always assume bounded potentials so that 
the conditions in Theorem 12.91 are satisfied. In the context of the density 
of states, we are interested in a Feynman-Kac formula for Hamiltonians on 
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bounded domains. Let us denote by A the set of all paths staying inside A 
up to time t, i.e., 

ft A = {b E C([0, t] : \b(s) G A for all < s < t} . 

Then e~ tH ^ is simply given by restricting the integration in 1)2. to the 
set Q A , e.g., 

A proof can be found in Simon |154j and Aizenman— Simon pQ. There is 
also a Feynman-Kac formula for Neumann boundary conditions (see |74j 
and references given there). 

Now we are able to state the probabilistic representation of the density 
of states measure in terms of Brownian motion. 

Theorem 2.10. The Laplace transforms of N(E), N D (E) and N N (E) 
agree and are given by 

N{t) = N D {t) = N N (t) = E x EW(e-ti Vu( X' y > da ). (2.12) 

For a detailed proof see, e.g., |68j and references there. The first step in 
the proof is to check that the right-hand side of 1)2. 12|) is finite for all t > 0. 
Interchanging the expectation values with respect to random potential and 
the Brownian motion, this follows from Jensen's inequality. 

The second step of the proof is to compare 1)2.12(1 and the Laplace trans- 
forms of the approximating density of states measure vl and . To prove 
the theorem one has to estimate the hitting probability of the boundary 
of A for a Brownian motion starting and ending far away from the bound- 
ary. Using standard facts of Brownian motion, this tends to in the limit 
|A| oo. 

Once we know that the Laplace transforms of N, N D and N N agree, it 
follows from the uniqueness of the Laplace transform that N, N D and N N 
agree themselves (see, e.g., |40j ). 

2.5. The Density of Surface States. We would like to define a den- 
sity of states measure for surface potentials as well. Suppose we have a 
surface potential of the form 

V*(x 1 ,x 2 ) = ^2 Qii(u)f( x i-h,X2) 

where, as above, x G M. d is written as x = (x\,X2) with xi G R ,X2 G M. d2 . 
In addition to the surface potential, there may be a random or periodic 
potential V b (x), which we call the "bulk" potential. The bulk potential 
should be stationary and ergodic with respect to shifts Tj parallel to the 
surface. Stationarity perpendicular to the surface is not required in the 
following. 

This allows "interfaces" in the following sense: Let d\ = d — 1, so the 
surface has codimension one. Thus it forms the interface between the upper 
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half space V+ = {x; X2 > 0} and the lower half space V-. The bulk potential 
V b may then be defined by V b (x) = V\(x) for X2 > and = for 
X2 < 0. Here V\ and V2 are random or periodic potentials on M. d . We set 
H b = H + V b , which we call the "bulk operator" and H u = H b + V£. 

We could try to define a density of states measure in the same way as 
in (|2.1j) . i.e., look at 

lim -Ltr(^(F w ) XA J. (2.13) 

L— >oo L 

It is not hard to see that this limit exists and equals 

E(tr( XAl ^(F b ) XAl )). (2.14) 

In other words, (|2.13|) gives the density of states measure for the bulk op- 
erator. After all, this is not really surprising. The normalization with the 
volume term L d is obviously destroying any influence of the surface poten- 
tial. 

So it sounds reasonable to choose a surface term like L dl as normalization 
and to consider 

lim J-tr(^) X Aj. (2.15) 

However, Definition (|2.15j) gives a finite result only when supp <pC\o(H b ) = 0. 

To define the density of surface states also inside the spectrum of the 
bulk operator, we therefore set 

v*{<P) = Hm -jjj tr ((<p(H u ) - ip(H b )) X A L ) • (2.16) 

Of course, it is not obvious at all that the limit (|2.16|) exists. In |36l 13 7| 
the authors proved that the limit exists for functions tp G Cq(IR). Hence 
the density of surface states is defined as a distribution. The order of this 
distribution is at most 3. Observe that, in contrast to the density of states, 
the limit in (|2.16j) is not necessarily positive for positive (p due to the sub- 
traction term. In fact, in the discrete case it is not hard to see that the total 
integral of the density of states, i.e., ^(1), is zero. Therefore, we cannot 
conclude that the density of surface states is a (positive) measure. 

Kostrykin and Schrader jl06L I107j proved that the density of surface 
states distribution is actually the derivative of a measurable locally inte- 
grable function. They do not prove that this function is of bounded varia- 
tion, thus leaving the possibility that u s is not given by a measure. See also 
the papers I17j by Chahrour for regularity properties of the density of 
surface states on the lattice. 

Outside the spectrum of H b , the distribution v s is positive, so that the 
density of surface states is a measure there. In jUj it was proven that below 
the "bulk" spectrum o~(H b ) the density of surface states can also be defined 
by using (Neumann or Dirichlet) boundary conditions. We expect this to 
be wrong inside the bulk spectrum. 
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3. Lifshitz Tails 

3.1. The Problem. For a periodic potential V the integrated density 
of states N(E) behaves near the bottom Eq of the spectrum o-{Hq + V) like 

N{E) ~ C(E-E ) d/2 . (3.1) 

This can be shown by explicit calculation for V = and was proved for 
general periodic potentials in |81j . 

On the basis of physical arguments Lifshitz \H7\ I118j predicted a com- 
pletely different behavior for disordered systems, namely, 

N(E) ~ d e ~ c ^ E - E ^~ d/2 (3.2) 

as E \ Eq > — oo. This behavior of N(E) is called Lifshitz behavior or Lif- 
shitz tails. The reason for this peculiar behavior is a collective phenomenon. 
To simplify the following heuristic argument, let us assume that V w > and 
Eq = 0. To find an eigenvalue smaller than E, the potential V w has to be 
small on a rather large region in space. In fact, to have an eigenvalue at 
small E > 0, the uncertainty principle (i.e., the kinetic energy) forces the 
potential to be smaller than Bona set whose volume is of the order E~ d / 2 . 
That is small on a large set is a typical "large deviations event" which 
is very rare — in fact, its probability is exponentially small in terms of the 
volume of the set, i.e., its probability is of the order 

e- c * E - d/2 (3.3) 

which is precisely the behavior (|3.2[) predicted by Lifshitz. It is the aim of 
this section to discuss the Lifshitz behavior (|3.2j) of the integrated density 
of states as well as its extensions and limitations. 

The first proof of Lifshitz behavior (for the Poisson model (|1.8|l ) was 
given by Donsker and Varadhan |33| . They estimated the Laplace trans- 
form N(t) for t — > oo using the Feynman-Kac representation on N (see 
Section 12. 4J) . Their estimate relied on an investigation of the "Wiener 
sausage" and the machinery of large deviations for Markov processes de- 
veloped by these authors. To obtain information about the behavior of 
N(E) for E \ = Eq from the large t behavior of N(t) one uses Tauberian 
theorems (HI II 1] . This technique was already used by Pastur |145l [B] and 
developed in |142l 11*4*] and recently in |125j . 

Donsker and Varadhan |33| needed in their proof of 1)3.2(1 that the single 
site potential / decays faster than (1 + |x|)~( d+2 ). They asked whether 
this condition is necessary for the result (Lifshitz tails) or just necessary for 
their proof. It was Pastur ( |146j ) who observed that, in fact, the Lifshitz 
asymptotic is qualitatively changed if / has long range tails, i.e., if f(x) ~ 
C(l + |x|)~ a for a < d + 2. Observe that a > d is necessary for the mere 
existence of V u . Pastur proved the behavior 
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as E \ Eq for d < a < d + 2. We call this behavior Pastur tails. For a 
disordered system with constant magnetic field in dimension d = 2, Pastur 
tails (JH2J were found for a// a > d = 2 in [10) . 

These results and more observations of the last several years indicate 
that the asymptotics of the integrated density of states even at the bottom 
of the spectrum is more complicated than expected. To be more precise 
and following the terminology of |148j , we can distinguish two qualitatively 
different behaviors in the low energy asymptotics of the integrated density 
of states. For short range potentials and "fat" single site distributions, 
the asymptotics of N(E) is determined by the quantum kinetic energy as 
predicted by Lifshitz. Hence it is called quantum asymptotics or quantum 
regime. On the other hand, for long range potentials or "thin" single site 
distributions, the leading asymptotics of the integrated density of states is 
determined by the potential, i.e., by classical effects. This situation is called 
the classical regime. 

We will discuss these phenomena in this section. We start with the short 
range case (quantum regime) . The proof of Lifshitz tails we present here is 
based on spectral theoretic arguments close to Lifshitz's original heuristics 
(see Section l3~2|) . 

We then discuss the long range case (classical regime) (Section I3.3|) 
to some extent, including recent results |86j of single site potentials with 
anisotropic decay resulting in a mixed classical-quantum regime fSection l3.4j) . 

Classical and quantum behavior of the integrated density of states and 
the transition between the two regimes is best understood for the Anderson 
model. The approach of [125] combines spectral theoretic and path integral 
methods. We will present this in Section fa. 51 

Lifshitz predicted the behavior (|3,lj) and IJM.2JI not only at the bottom 
of the spectrum but also for any band edge of the spectrum. To distinguish 
these two cases, we will speak of internal Lifshitz tails in the latter case. 
Investigating Lifshitz behavior at internal band edges turns out to be much 
more complicated than at the bottom of the spectrum. In fact, already 
the investigation of periodic potentials at internal band edges is extremely 
complicated. We will discuss internal Lifshitz tails (following |97H104|ll05j ) 
in Section l3~THl 

Finally, we will look at random Schrodinger operators with magnetic 
fields in Section EH1 



3.2. Lifshitz Tails: Quantum Case. 

3.2.1. Statement of the main result. The aim of this subsection is a proof 
of Lifshitz behavior close to his original heuristics and without heavy ma- 
chinery. We will prove the quantum asymptotics in l|3,2j) for short range 
single site potentials and "fat" single site distributions. We will make no 
attempt to reach high generality but rather emphasize the strategy of the 
proof. 
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As before, we consider random alloy-type potentials of the form 

VUx) = J2 Qi(u)f(x-i). (3.5) 



We assume that the random variables qi are independent and identically 
distributed with a common probability distribution Pq. We suppose that 
the support of Po is compact and contains at least two points. 

As always, we also suppose that the single site potential / is non- 
negative, bounded and decays at infinity as fast as | x\~( d+£ \ The technique 
we are going to present allows us to treat local singularities of /. (We refer 
to HOI US] for details.) 

To ensure Lifshitz tails in the sense of 1)3.2)1 . we need two conditions: 

Assumption 1: Define q m %n = inf supp(Po). We assume that 

Po([ Qmim 9min + e)) > Ce N (3.6) 

for some C, iV and all e > small. 

Condition (jHSJ) means that the distribution Po is "fat" at the bottom of 
its support. Note that this condition is, in particular, satisfied if Po has an 
atom at q m im i- e -i if Po({feun}) >0. 

The second condition we need is precisely the "short range" condition 
already encountered by Donsker and Varadhan |33j . 

Assumption 2: 

f{x)<C{l + \x\)" {d+2 \ (3.7) 
We are ready to formulate the main result of this subsection. 

Theorem 3.1. If Assumptions (1) and (2) are satisfied, we have 

ln(-lnMP)) d , , 

lim — K -— \-4 L = — • (3.8 

e\e ln(P-Po) 2 v ' 

Observe that equation (|3.8ft is a weak form of Lifshitz's original conjec- 
ture (|3.2j) . In their work |33j . Donsker and Varadhan proved the stronger 
form for the Poisson potential 

lim ^ = -C d (3.9) 

where Cd is a (computable) positive constant. 

Both the short range condition (Assumption 2) and the fatness condition 
(Assumption 1) turn out to be necessary for the above result, as we will see 
later. For example, if the single site potential / decays substantially slower 
than required in the short range condition the integrated density of states 
decays faster than in the 1)3.8)1 . 

We define the Lifshitz exponent 7 by 

ln(-lnMP)) 
7= lim ^ 3.10 

e\e ln(E-E ) 
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whenever this limit exists. With this notation we may rephrase (|3.8jl as 
7 = — d/2. The Lifshitz exponent for periodic potentials is 0. 

3.2.2. Strategy of the proof. The proof of Theorem 13.11 consists of an 
upper and a lower bound. The next subsection will provide us with the 
tools we need for these bounds. 

It will turn out that the bounds are easier and more natural for posi- 
tive random potentials. Therefore, we will split the random potential in a 
periodic and a positive random part 

V u (x) = ^ <?min f{x - i) + ^ (lii^) ~ Imin) f(x ~ l) 

= V pcr (x) + VUx). (3.11) 

We will subsume the periodic potential under the kinetic energy and denote 
the positive random potential V w in a slight abuse of notation again by V^. 
Thus we have 

H w = H l + V u (3.12) 

with Hi = Hq + V per a Hamiltonian with periodic "background" potential 
Vp er and 

where the independent qi > have a common probability distribution Pq 
with = inf (suppPo)- 

For the upper bound below, we need information about the two lowest 
eigenvalues of Hi restricted to a box. If V pcr = 0, these eigenvalues can 
be computed explicitly. However, if V pcT ^ 0, we need a careful analysis 
of periodic operators. This was done in |81| and |129j . Here we restrict 
ourselves to the case V pei = 0, avoiding some technical complications. Note 
that this implies Eq = inf (a(H u )) = 0. We refer the reader to the papers 
|80j and |129j for the general case. We also remark that |86j contains an 
extension of the approach presented here that works for Poisson potentials 
and various other potentials as well. 

3.2.3. The Dirichlet-Neumann bracketing. The first step in the proof is 
to bound the integrated density of states from above and from below using 
the Dirichlet-Neumann bracketing as in Corollary 12.71 We have 

-L. E {N(H° Al ,E))< N(E) <j^jE {N(H^ Al , E)) . (3.14) 
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The side length L of the cube Al will be chosen later in an ^-dependent 
way when we send E to Eq. We estimate the right-hand side of 1)3.1 4 J) by 

E(N(H^ Al ,E) = j N(H^ Al ,E) dF 

= I N(H^ Al ,E) dF+ ( N(H^ Al ,E) dF 

Je 1 (h^ Al )<e Je 1 {h* Al )>e 

<F(E 1 (H^ Al )<E) N(H ^ l ,E). 

With N(H ^ L ,E) < (Ci + C 2 E) d / 2 \A\ following from Weyl asymptotics, 
we get for < E < 1 the estimate 

-L P (^(^aJ < E) < N(E) <CF {E 1 (H^ Al ) < E). (3.15) 

The problem now is to find upper and lower bounds such that after taking 
the double logarithm, the left- and the right-hand side of 1)3.15)1 coincide 
asymptotically. In general the upper bounds are more difficult than the 
lower bounds. To prove the lower bound we only have to "guess" a good 
test function, whereas for the upper bounds, one has to prove that all eigen- 
functions for energies in [0, E) roughly behave the same way. 

It is an astonishing fact that the lower bound from 1)3.15)) in all known 
cases leads to the asymptotically correct behavior of the integrated density 
of states, a fact emphasized by Pastur. 

3.2.4. The lower bound. For simplicity we restrict ourselves to single site 
potentials / with supp / C Ai so that /(• — i) and /(• — j) do not overlap for 

i 7^ j. The necessary changes for the general case will become clear when 
we discuss long range potentials /. 

By the Neumann-Dirichlet bracketing in 1)3. 15|) . we have for arbitrary L 
and any tp G D(Aj>J with \ \^\\ L ^{A L ) = h 

N(E) > \Al\- 1 !? (EtiH^J < E) 

> \Al\- 1 F((^,H^)l2 {Al) <E). (3.16) 

A natural choice of ip for 1)3.16)1 seems to be the ground state 'ipo of —A A , 

ipo(x) = |Ax| — 2 . Unfortunately, this function does not obey Dirichlet 
boundary conditions and is therefore not admissible for 1)3.16)1 . 

This problem can be circumvented by multiplying tpo by a function which 
is zero at the boundary of Ao- To do so, let us take x £ C°°(W i ), suppx C Ai, 
x{ x ) = 1 on Ai and < x( x ) — 1- We set xl(x) = x{j) an d = 

Xl(x)Mx)- Then ^ L G D(Aj> L ), \\^ L \ \ > \ and 
(tl>L,H u il> L ) < (4> ,H w i; } + CL- 2 

= \Al\~ 1 [ V^dx + CL- 2 . 
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Note that the "error term" L 2 is due to the influence of the kinetic energy 
(a second-order differential operator). Inserting in (|3.16j) we get 

N{E) > \A L \~ l P (\A L \- 1 jf V u [x) <E- Or*} 

> \Al\~ 1 P^IAlI" 1 ii / ||i Y <&M < E ~ CL ~ 2 ) • ( 3 - 17 ) 

In principle, we can choose L as we like. However, if E < CL~ 2 estimate 
(|3.17j) becomes useless. So it seems reasonable to choose L = (3E~z and we 
obtain 

(E33) > IAzT 1 p(Wl _1 Y ftM < ce 

> \Al\~ 1 P (g < C E) Ld 

> dE-i (C E n )° 2eA . 



Thus we conclude 



M-\n(N(E))) d 
lim ^ - V ,!> > --. 3.18 

E\o InE 2 y ' 



3.2.5. The upper bound. The strategy to prove the upper bounds of 
P (E\{H~f A ) < E) in Q3.15[) can be divided in two parts. The first step 

is to find a lower bound for Ei(H® A ) in such a way that it is possible to 
control the influence of the random potential. This is done by an application 
of Temple's inequality. The second step is to balance between the size of Al 
and the probability of a random potential such that most of the potential 
values are small. 

We start by stating Temple's inequality for the reader's convenience. A 
proof can be found, e.g., in |150j . 

Theorem 3.2 (Temple's inequality). Suppose H is a self-adjoint oper- 
ator, bounded below which has discrete spectrum and denote by E n (H),n = 
1,2,... its eigenvalues (in increasing order, counting multiplicity). If fx < 
E2(H) and ip € D (H) with \\ip\\ = 1 satisfying {tp,Hip) < fj,, then 

E l iH ) >(>i>,m,)- { < p ' H ' ! t :) - {i '' H * ? 



To apply Temple's inequality, we set £b(— ^a l ) := l 1 — ^{H^^ ). Note 
that by direct computation, Ei(—A% ) = = E and E 2 {-A% ) ~ L~ 2 . 
These facts require a careful analysis if there is a periodic background po- 
tential as in (JTTTj) and (|3~T2j) : see [HT|. 

Next we need a good approximation ip of the ground state of Huj\ l • This 
is done by choosing to be the ground state of — Aj[ , which is intuitively 
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close to the correct ground state for small E. The function ijj is given by 
?P(x) = IAlI" 1 / 2 . 

To apply Temple's inequality we have to ensure that with the above 
choice, (ip,Hijj) < /i ~ cL~ 2 . We force this to happen by changing the 
coupling constants qi to qi = min (qi(u>), aL~ 2 ) with a suitable a > 0, small 
enough. If H denotes the corresponding operator, we have E\(H) > E\(H). 
An application of Temple's inequality to H and an elementary calculation 
yield the following lemma. 

Lemma 3.3. 

Exilic > I -L £ (3.19) 

«eA L 

A consequence of the lemma above is the intuitively convincing estimate 
F(E 1 (H U ^ L )<E) <p(|^-t E 9i(")^ 2E ) ( 3 - 20 ) 

^ P (j7 E 9iM<2£j. (3.21) 

The expression ()3.21(l for E small very much resembles a large deviation 
probability which would lead to a bound exponentially small in the volume 
term L d . At first sight, Cramer's theorem, a result of the theory of large 
deviation, seems to be applicable (see, e.g., |3H I32|. I58j ). 

However, there is a complication here: To obtain a large deviation event 
in (I3~2T1) we need that E(&) < E. Thus, if we set L = L(E) = /1ET 1 / 2 
with P > small, the event (|3.21|) is, indeed, a large deviation event and we 
obtain the following bound. 



Lemma 3.4. 

TP I 

IA, 



Ll i€A L ' 

for E close enough to zero and L < (3 E~ l l 2 . 

Combining the results above, (|3.15j) and (|3.2()j) . we have proven 

N{E)<C^ E ~ d '\ 

3.2.6. Final remarks. The idea of using Neumann-Dirichlet bracketing 
to prove Lifshitz tails first appeared in |77| . It was carried over to the dis- 
crete Anderson model by Simon |156j who streamlined it at the same time. 
The proof was extended to more general alloy-type potentials by Kirsch and 
Simon |80j . who still needed reflection symmetry of /. Mezincescu |129j 
modified the upper bound by introducing other boundary conditions to get 
rid of this extra assumption. We refer to |86j for a rather general proof 
using these techniques. 
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3.3. Long Range Single Site Potentials: A "Classical" Case. In 

this section we turn to an example of classical behavior of the integrated 
density of states near = inf(cr(-£^)), in the sense of Section 13.11 namely, 
to long range single site potential /. 

The upper bound on N(E) is easier than for the short range case. While 
there is a subtle interplay between the kinetic energy and the potential in 
the short range case (/(x) < |x|~^ +2 ^), it is the potential energy alone that 
determines the leading behavior of N(E) (E \ 0) in the long range case. 

Assumption: In this section we suppose that 

< /(*) < „ , C , nn . (3.22) 



{l + \x\) a ~ J w - {l + \x\Y ' 
Theorem 3.5. Assume (|3.fi|) . // (|3.22j) holds for an a with d < a < 



d + 2, then 



e\e InE a - d 



In the terminology of (|3.1U|) . Theorem 13.51 states that the Lifshitz expo- 
nent for the long range case (a < d + 2) is d/(a — d). 

Proof. To simplify the argument, we assume as in the short range 
case, there is no periodic background potential and q m i n = 0. Consequently, 
Eq = 0. We start with the upper bound and estimate 



Bk(H u l)> inf £ 



zeAi (x + |i| 



ie. 



Hence 



N(E)<C 1 F(E 1 (H^ 1 )<E) 

<c 1 f(y qi - < e] 



<Ci 



\i\<L 

^ < C 3 EL^A . (3.24) 

\i\<L J 



8 l 

We choose L = -^rE a - d {(3 small). Hence 

V |<|<£ 

If f3 < 2 E(go), standard large deviation theory gives 



I I T ' 



CL d _ P -CE 

L d 

\i\<L 
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We turn to the lower bound. As in the proof of the lower bound QH.17JI 
in the previous section, 

N(E) > -r^P f j^-r [ V u (x)dac <E- CL' 2 ] 
\^l\ Ja l J 

^ih P (Y^TTl I f(x-i)dx<E-CL- 2 ). 
|Az| \f^ d | Ai| Ja L ) 

Due to the long range tails of /, we cannot ignore the summands with \i\ 
large. Instead, we estimate 

Y * TT~r / f( x ~ l ) dx 

- TTT Y * f(y) d V+ Y 2*Ta~T / f(x-i)dx 

|ALl |i|oc<2L J |i| 00 >2L |AL|J ^ 

- T/T~~\ Y + ^ max X] Ta~7 / ~ *) 

|A2L| |i|Sui M ^2I |Ai| ^ 

where g max = sup(suppPo)> Pq being the distribution of the We estimate 



\i\>2L 1 1 JAl \i\>2L 1 1 17 Al 1 

|i|>L 

<c 6 



-dx 



L a-d- 



Thus we obtain 



N(E) > T^ P (]^7 Y * ^ °7 E ~ C » L ~ 2 ~ W a ~A • (3-25) 

i|<2L 

As the derivation shows, the L~ 2 term comes from manipulating the kinetic 
energy, while the L~( a ~ d ^ term is due to the potential energy. Note that for 
a > d + 2, the term L~ 2 (kinetic energy contribution) dominates in (|3.25|) . 
In this case we can therefore redo the estimates of Section T3.2.4I and obtain 
a lower bound as we got there. However, for a < d + 2, the term L~( a ~ d ^ 
wins out in (|3.25|) . Remember, this term is due to the potential energy 
distribution. We obtain 

N(E) > -Lpf-L- V q, t < C 7 E - C 10 L^ a ~ d A . 
\A L \ V|A2l|^ l J 
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This time, E has to be bigger than L~( a ~ d ^; more precisely, E > Ci\L~^ a ~ d \ 
i 

Hence L = C\ 2 E a ~ d so 

N(E) > ( qi = for \i\ < 2L) 

I Al| 

" |Al| c 



d — 



> C u E^i e~ Cl5E . □ 

3.4. Anisotropic Single Site Potentials. Recently, Theorems 13.11 
and !3.5l were generalized to single site potentials / decaying in an anisotropic 
way at infinity (jSp). Let us write x G R d = R dl x P^ 2 as a; = (x 1 ,x 2 ), 
x\ G M^ 1 , X2 G M d2 and suppose that 



i — i i — r- < fix) < i — i i — i— ( 3 - 26 ) 

+ \x 2 \ a2 ~ \xi\ a i + \x 2 \ a2 

for \xi\, |a?2 1 > 1, and define in the usual way 

= £>/(*-*)■ (3-27) 

Let us define 7j = — and 7 = 7i + 72- Then the sum in 1)3.27(1 converges 
(absolutely) if 7 < 1. 

In |86j the authors prove that there is Lifshitz behavior of N(E) for 
potential as in (|3.26|) and (|3.27|) in the sense that the Lifshitz exponent 77, 
defined by 

ln\ln(N(E))\ 

77 = hm — 

E\E In E 

exists (Eq = inf a(H ul )). r/ depends on the exponents Ckj, of course. If both 

J»_<* an d _^<^ (3.28) 
l-7~2 I-7-2 v ; 

we obtain the "quantum" exponent: 

d 

Observe that (|3.28|) reduce to the condition a < d + 2 for the isotropic 
case «i = Q2 = a. If 

71 >^ and ^ > d l (3.29) 



1-7 2 I-7 2 

we are in the "classical" case both in the d\- and ^-directions. Then 

7 

V = ~z • 

1 - 7 

The third case 

71 . di 72 ^ d 



<^ and -^>^ (3.30) 



I-7-2 I-7 2 
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is new compared to the isotropic case. It is, in a sense, a mixed quantum- 
classical case. The Lifshitz exponent is given by 

d\ 72 



2 1-7 

We note that the cii-direction and the ^-direction "influence each other" in 
a rather sophisticated way. In |86j these results are proved for alloy-type 
potentials as well as for Poisson (and related) models. We summarize: 

Theorem 3.6 (HSJ). Suppose (|3~H)) and (|3~^)) hold. Set ^ = di/on and 
7 = 71 + 72 . Then the Lifshitz exponent n is given by 

J4 71 I f d 2 72 \ , q on 

r? = — max < — , > — max < — , >. (3.31) 

' I 2 ' 1- 7 J I 2 ' 1 - T J V j 

3.5. Path Integral Methods and the Transition Between Quan- 
tum and Classical Regime with Respect to the Single Site Mea- 
sure. We start this section with two observations which indicate that the 
asymptotic behavior of the integrated density of states depends qualitatively 
on the distribution Pq of the eft. 

Let us first assume that the "fatness" condition (|3.ti|) is satisfied in the 
strongest form, namely, inf supp(Po) = and P(qo = 0) = a > 0. An 
inspection of the proofs in Section 13.21 shows that for this case we have 
actually proven 

lim inf lnN ( E \ /o > -d, (3.32) 
(E-E )- d / 2 ~ y J 

In N(E) 

UP (E-E )-W ~ 2 (3 " 33) 

with d,C 2 > 0. If instead P(q = 0) = (but still P(q > e) > Be n ) the 
lower bound requires a logarithmic correction 

liminf ^ N f] > - Cl . (3.34) 

{E - Eo)-*/*\]n(E - Eo)\ ~ V ' 



The second observation concerns unbounded single site measures Po- bi 
it is proved that both the classical and the quantum regime can occur for 
the discrete, unbounded Anderson model and more general matrix opera- 
tors. Depending on the single site measure, collective phenomena may occur 
similar to those we encountered above. In other situations, the single site 
measure alone determines the behavior of the integrated density of states. 

It seems difficult to understand the mechanisms causing the transition 
from quantum to classical regime with respect to the single site measure by 
using the spectral analytic approach close to Lifshitz's original intuition. 

The first approach to prove Lifshitz tails is based on the Donsker- 
Varadhan technique (see also Section 13.11) . This method to compute the 
Laplace transform of the integrated density of states in the limit t — > 00 
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is a far-reaching generalization of the Laplace method known from classical 
analysis. The starting point is the path integral representation 



The Donsker-Varadhan technique is based on a large deviation principle 
satisfied by the product probability measure dP x cflPg q combining the ran- 
dom potential and the Brownian motion. In an informal sense, it makes it 
possible to quantize the asymptotic probability of a Brownian particle to 
stay most of its lifetime in a pocket with a favorable configuration of po- 
tential values. Using the large deviation principle, one can balance between 
favorable configurations and their small probability by applying Varadhan's 
lemma. Last but not least, given the large time asymptotics of the Laplace 
transform N(t), one can reconstruct the Lifshitz tail behavior using Tauber 
theory. 

The Donsker-Varadhan technique was worked out by Nakao |142j for 
the Poisson model with / > 0. He proved 

InN(E) 
hm J7-- = — Cd 

E\0 E~ d / 2 

where Cd is a (computable) positive constant. In the 1990's Lifshitz asymp- 
totics became a starting point for stochastic analysis of diffusion in random 
media. We mention the work of Sznitman (see, e.g, |167^ 11681 11691 H7()j ) 
in the continuous case, especially for Poisson potentials, and in the discrete 
context the moment analysis for the so-called parabolic Anderson model 
(PAM) (see, e.g., [7| ,I48^H9^I5()) ). Here Brownian motion has to be replaced 
by the continuous time Markov chain generated by the discrete Laplacian. 

The phenomenology described at the beginning of this subsection was 
also observed in the moment analysis of the parabolic Anderson model start- 
ing in |48j . As we will see, the case of the double exponential distribution 
discussed in |49j can be interpreted as the borderline between the quantum 
and the classical regime. The paper clarified the discrepancy between 
the lower and the upper bounds in 1)3.32(1 , 1)3.33(1 and 1)3.34(1 . Still, a general 
principle explaining the transition from quantum to classical regime with 
respect to the single site measure was not formulated. 

We want to systemize the phenomenology discussed above in the follow- 
ing theorem taken from |125j . To combine the bounded and the unbounded 
case, we assume that the cumulant generating function is finite, e.g., 

G(t) := log E (exp(-tK,(0))) < oo (3.35) 

for all t > 0. Furthermore, we apply the Legendre transformation to define 
the rate function 

1(E) := sup[Et - G(t)] (3.36) 

and for t, A > 0, we set 

S(X, t) := (Ai)~ 1 G(At) - t _1 G(t). (3.37) 
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Informally the scale function S(X, t) measures the change of the cumulant 
generating function after rescaling the time. 

Theorem 3.7. We consider the discrete Anderson tight binding operator 
hu> = hq + Vu and set Eq = inf a(hu,). Suppose G(t) < oo for all t > 0. Then 
we distinguish the following four cases: 

(i) Let S(X,t) ~ c (X p — l)t p with c,p>0. Then the IDS behaves in the 
limit E — > Eq = — oo like 

log(N(E)) = -I(E + 2d + o(l))(l + o(l)). 

(ii) In the case S(X, t) ~ c log(A) with p = and c > 0, we have in the limit 
E — > Eq = —oo, 

-KI(E + C 2 )(l + o(l)) < log(JV(^)) < + Ci)(l + o(l)) 

with K > 0, 

and 

a(d) 



C 2 = if(rf) min 



-c + clog(c 1 ),max 



1 - 4exp(-iTc), 



(hi) In the case S(X,t) ~ c(l - A^)t^ tuifft -1 < p< and c> 0, t/ie ZDS' 
behaves in the limit E \ Eq = like 

- J ftTi£;- 1 /2(d-2p- 1 (p+i))( 1 + ( 1 )- ) < io g (jv(£;)) 

< -Jf 2 i?™ 1 / 2 ( d - 2 ^ 1 ^+ 1 ))(l + (l)). 

(iv) In the case S(X,t) ~ — c (At) -1 log(i) wii/i c > 0, we have in the limit 
E\E = 0, 

—K\E~ d l' 1 \°&{E)(\ + o(l)) < log(N(E)) 

< -if 2 S- d / 2 log(E)(l + (l)). 

The scaling assumption S(X, t) ~ c (A p — l)t p with c, p > in the first 
case corresponds to "fat" unbounded single site distributions; the behavior 
of the integrated density of states is classical. The second case represents the 
double exponential case, while in the third situation the single site distri- 
bution is bounded, but very thin. The fourth case corresponds to relatively 
fat single site distributions studied in |74j . 

Although, by now, there are results covering a lot of possible single site 
distributions, there seems to be no systematic approach known to explain 
this phenomenology. Furthermore there exists two relatively different ap- 
proaches as discussed above. A first step to combine the functional analytic 
and the path integral approach as well as to systemize the known results 
with respect to the single site distribution seem to be [T251IT26] . In contrast 
to the direct analysis of the operator H u in Section 13.21 but in analogy to 
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the path integral methods, one is interested in the large time behavior of the 
semigroup exp(— h u t). The Lifshitz asymptotics follows by an application 
of (modified) Tauber theorems |127j . 

The first step of the argument in |125j is to restrict to a (time- 
dependent) box A = A t (0) by introducing discrete Dirichlet boundary con- 
ditions and to approximate E[exp(— h^t)^, 0)] in the limit t — > oo by 



E[exp(-^(^) t )] = E 



sup exp (-t [{y/p\h% VP) + (Vp\V(w)y/p)]) 
.peMi(A) 

(3.38) 

Here Ei(h^) = mfa(h^) is the principal eigenvalue of and Mi (A) is 
the set of probability measures on A. Equation (|3,38|) is a consequence of 
the min-max principle and the nonnegativity of the ground state. It is the 
starting point to find upper and lower bounds of the Laplace transform of 
the IDS. 

To illustrate the central effect explaining the transition from the quan- 
tum mechanical to the classical regime, we want to sketch the very elemen- 
tary proof of the lower bounds starting from (|3.38|) . The first step is to 
interchange the expectation value and the supremum 

E [eM-Ei(h^)t)] = E sup exp (-t [{y/p\h%y/p) + (^l^)v^)]) 

LpeMx(A) 

> sup exp(-t (v / p|^ov / P)) E ex P( "^PW^t 15 )* 



pGMi(A) 



= sup exp l-t (VPI/4VP) +5^G(p(x)t) 
peAh(A) V zeA 

(3.39) 

The second step is to define a subset D C Mi (A) of relatively uniform 
probability distributions concentrated on a subvolume of A. With the side 
length L of A, we set 1 < I < L and A; := {x € Z d : \x\oo < I}. The ground 
state of the discrete Laplacian h 1 restricted to A; with Dirichlet boundary 
conditions is given by 

{ : Aj-> [0,oo), (3.40) 

♦-n(^r)" !s KS) < 3 - 4i > 

j=i 

and the corresponding principal eigenvalue of h^ 1 is 

E l (h^) = 2d(l-co B (-^\\ (3.42) 

The subset D C Mi (A) is then defined by 

D := {$ : 1 < I < L}. (3.43) 



IDS FOR RANDOM SCHRODINGER OPERATORS 



29 



Restricting the estimate to D, we get 

E(exp(-t E 1 (h*))) 

exp(G(t)+t sup ( -(^\h^) + Y,P(^)S(p(x),t) 



> 



peMi(A) 



exp(G(t)+t supf -(VpKVp) + y2p{x)S{p(x),t) 

Using the definition of S(\,t), the convexity of the cumulant generating 
function G(t) and the Jensen inequality, we can estimate for p G Mi(Aj) C 
Mi (A), 



< £ G(p(x)t)) - r 1 ^*) 



xeA x&Ai 

= rH d [ , 



> t-H d G(r d t p(x)) - t -1 G(t) 

V _f A. / 



XGAj 

_ G(r d t) _ G(i) 
= S(r d ,t). 

So the uniform distribution on A; minimizes ^2 xe \P(x)S(p(x), t) with re- 
spect to the variation over p G Mi(Aj). We have 

E(exp(-t£i(#f(w)))) >exp(G(t)+t sup f-4dsin 2 f J — L V 5 ^, *) 

V 1<1<L V ^ ' + 1/ 

and the only remaining problem is to maximize 



-4^sin 2 (| J--) + i<?(r d ,t) 



' 27TT 1 



' c{l~ d P -l)tP +1 p>0 
-cdt log(Z) p = 

c (1 - l~ d P)tP +l -l<p<0 
k -cZ d log(i) 

with respect to Z. The exponent of the time i in the scaling expression 
is responsible for the occurrence of the classical or the quantum regime. 
In the case p > 0, the scaling term increases faster in t than the linear 
time dependence in the diffusion term. Consequently, the maximum will be 
asymptotically 1 = 1. This corresponds to the classical regime. In the case 
— 1 < p < as well as in the fourth situation, the scaling term is sublinear 
and the diffusion term is dominating. Like in Section I3.2| a collective be- 
havior of potential values is necessary and we are in the quantum regime. 
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In the case p = 0, the diffusion and the scaling term are both linear in time. 
So the optimal peak size depends strongly on the constants. This is the 
borderline between the classical and the quantum regime. It corresponds to 
the double exponential distribution. 

The upper bounds are much more complicated. It is not possible to 
interchange the supremum and the expectation with respect to the random 
potential. Moreover, we have to to estimate 1)3.38(1 for all p G Mi (A). The 
first problem is solvable by a variant of the ordinary Laplace method. The 
second problem is attacked using the convexity of the cumulant generating 
function G(t) and ideas from spectral geometry. For details, we refer to 
[HU and P2H]. 

3.6. Internal Band Edges. Lifshitz predicted the "Lifshitz behavior" 
not only for the bottom of the spectrum but also for other band edges. 
We refer to this phenomenon as "internal Lifshitz tails." Internal Lifshitz 
tails have been proven for the Anderson model by Mezincescu |128j and 
Simon |157| . Their proofs apparently cannot be translated to the continuum 
case. In fact, the band edges of the Anderson model which they can handle 
are those coming from gaps in suppPo together with the boundedness of 
the kinetic energy. (To be more precise: Since for the Anderson model, 
II ^o|| 4c£ and /io > 0, there are gaps in the spectrum whenever there are 
gaps in supp Po of length exceeding Ad) . 

One can also handle the case of a point interaction potential in one di- 
mension, a problem which essentially reduces to a lattice problem. Formally 
this potential is given by 



where 5 is the Dirac- "function." This potential is also known as the random 
Kronig -Penney model. It turns out in this case that the lower edges do and 
the upper edges (for qi > 0) do not show Lifshitz behavior but polynomial 
behavior of A?" as for periodic potentials |79j . This is due to the fact that 
the upper edges are "stable boundaries" in the sense of |148j . The case of 
general one-dimensional alloy-type potentials was treated in |13()j . 

The multidimensional case is by far more difficult. The reason is mainly 
that periodic potentials are much less well understood in higher dimensions. 
For example, it is not true in general that bands are parabolic, as is the case 
for d = 1 and for the ground state band in arbitrary dimension. 

The paper |97j marks a breakthrough in this topic. Klopp uses the 
method of approximation by periodic potentials. Compared to Dirichlet- 
Neumann bracketing, one loses monotonicity, a property which was very 
useful above. However, Klopp manages to prove an exponential convergence 
rate for the periodic approximations. 

As mentioned above, not so much is known about the behavior of the 
band functions (of the periodic operators) at internal band edges. In fact, 
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Klopp has to make assumptions on the behavior of the integrated density 
of states for the periodic operator. 

Like Lifshitz tails at the bottom of the spectrum internal Lifshitz tails 
can be used as an input for a localization proof |175j . 

We consider an alloy-type potential with a continuous single site poten- 
tial / > 0, not identically equal to 0, with decay 

f(x) <C{l + \x\)-^ d+2+e \ (3.44) 

The random coupling constants are independent and have a common proba- 
bility distribution Pq with q m \ n = inf suppPo- We set V peT = ^iez d ^min f(x— 
i) and denote the integrated density of states of H pcr = Ho + V pcT by 
N peT (E). Furthermore, we suppose that E- is a lower band edge of H pcv , 
i.e., E_ E cr(H pcr ), but (£L - a,EJ) n a{H per ) = 0. 

It is reasonable to assume that generically iVp er behaves like (E — E-) d / 2 , 
for E \ E-, as it would for a unique parabolic band. In fact, this behavior 
is known in one dimension and for the bottom of the spectrum in arbitrary 
dimension. However, it is not clear that this is true in general, even not 
generically (see, however, |l()3j ). 

Thus, we have to assume such a behavior of iVp er : 
Assumption: Suppose i?_ is a lower band edge of H pcT . We assume that 

E\E- ln(E-E.) " ~ 2 ' {SAb) 

Under this assumption, the main result of |97j is: 

Theorem 3.8 (Klopp). If assumption ()3.45j) holds, then 

In (-In (N(E)-N(E_))) d 
lim — ^ J - = --. (3.46) 

For the case d = 2, one has more information about the periodic opera- 
tors |105j . In particular, there is always exponential decay of the integrated 
density of states at band edges. We refer to the review |100j for an intro- 
duction and further results. 

3.7. Lifshitz Tails for Surface Potentials. In this section we con- 
sider surface potentials of the form 

V£(x 1 ,x 2 )= ^2 qii{u)f(xi - h,x 2 ) (3.47) 

and suppose we have some spectrum below 0. This is the case if q m \ a = 
inf suppPo is negative enough. Note that for d 2 < 2, there is negative 
spectrum as soon as q m \ n is negative. For d 2 > 3, there is a threshold 
7 > such that the spectrum starts at if q m i n > —7 and there is negative 
spectrum if q m \ a < —7. 

We are going to investigate Lifshitz tails for surface potentials with Eq < 
0. Below the bulk spectrum (which starts at 0), the density of surface states 
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is positive, hence a measure. We may therefore define the integrated density 
of surface states N S (E) to be the corresponding distribution function. 

As before, we decompose the potential into a non-random background 
potential and a positive random potential 

V*(xi,x 2 )= ^2 qi 1 (uj)f(x 1 -i 1 ,x 2 ) 



iiez d i 



grain /(^l - h,x 2 ) + ^2 (Qh - Qmin) f(xi - ii,x 2 ) 



= V° p (x) + FJ(x). (3.48) 

The Neumann-Dirichlet bracketing technique goes through for this case 
as soon as we have sufficient knowledge about the background operator 
Hi = Hq + V s s p . Since we want Eq < — which makes the bottom of the 
spectrum "surface spectrum" — there is no case V sp = here. Moreover, the 
background potential is only periodic for the di-directions, but decays 
perpendicular to them. 

The analysis of those partially periodic potentials and the Lifshitz esti- 
mates for surface potentials were done in |87| for the continuous case. We 
assume that 

Po{[ gminj gmin + £)) > Ce N 

and 

< f(x u x 2 ) < f (l + \xi\)-^ +2 \ 
We also assume that f(xi,x 2 ) decays uniformly in ^-directions. Then we 
have: 

Theorem 3.9. // E < 0, then 

M -MN.W)) = _,h, 

e\e ]n(E-E ) 2 v ; 

There is also an analogous theorem for long range /. The paper |72j 
proves Lifshitz tails for surface potentials in the discrete setting by fairly 
different techniques. This paper also contains an analysis at the energy 
E = 0, i.e., for surface corrections to the bulk Lifshitz tails. 

3.8. Lifshitz Tails for Random Landau Hamiltonians. We turn 
to the density of states for operators of the form 

H w = H (B) + V u 

with a constant magnetic field B > and a non-negative random potential 
V w . (For a careful definition of the density of states and some basics, see 

[US ECU QUI ) 

We discuss the two-dimensional case first. The Landau Hamiltonian 
Hq{B) is given by 
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Hq(B) has a pure point spectrum for B ^ and d = 2. In fact, the 
eigenvalues are given by the "Landau levels" (2n + 1)B; n S N and all 
Landau levels are infinitely degenerate. One possible choice of the ground 
state is 

fofx) = - e -f |x|2 (3.50) 

7T 

which will play a major role below. For V w we take a Poisson potential or 
an alloy-type potential with q m \ n = 0. In this case, the bottom Eq of the 
spectrum of H w is given by the lowest Landau level, which is B. In jlOj the 
authors proved for the Poisson model the following result. 

Theorem 3.10. If B / and 

the Lifshitz exponent rj for H (B) + V W is given by 

V=— (=^-) (3-52) 
' 2-a V d-aj v ; 

for all a > d. 

This means that, according to our classification above, we are always in 
the classical case for d = 2 and constant magnetic field. 

Proof. We sketch the lower bound only and restrict ourselves to the 
alloy-type case. As usual we have to estimate 

P(E 1 (H UJ ^ L ) < E + E) 

from below. This time we have a ground state tpo for Hq(B) which is I? . We 
modify ipo near the boundary of Al to make it satisfy Dirichlet boundary 
conditions. Due to the (super-)exponential decay of ipo, the error we make 
is of the order (at most) e~ Co . Thus 

P(£i (H^J <E + E) >FN V u (x)\^\ 2 dx<E-e- CoL ' 

>P^y V UJ (x)dx <dE- e - C2L2 
= (*) • 

At this point we can literally repeat the estimates in the proof of Theorem 
13.51 and obtain in analogy to (|3.25j) 

(*) > P(TT^-r J2 Qi< CiE - C 2 e~ cL2 - C 3 L~^) . (3.53) 

|A2l| 

The only difference to the previous case is the error term due to the kinetic 
energy. In Theorem l3.5l it was of the order L~ 2 causing the different behavior 
for a > d + 2 and a < d. In (|3.53|) the error term is exponentially small, 
thus being negligible with respect to the potential term L^ "^ for all a. 
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Consequently, we may choose E ~ L a - d . By a large deviation estimate, we 
obtain 

(13331) > P(ft < C'E)^ 
> Me~ dLd 

= Me- dE ~^~ d . □ 

Theorem I3.1UI implies that for compactly supported /, the integrated 
density of states N(E) decays subexponentially. 

Erdos |38l I39j proved that it decays, in fact, polynomially. Erdos' proof 
is based on a careful estimate of the Laplace transform of N. It uses an 
analog of the Feynman-Kac formula for magnetic Schrodinger operators, the 
Feynman-Kac-Ito formula (see |154j ) . There are a couple of complications 
in the Feynman-Kac expression of N due to the magnetic field. The most 
serious one is the fact that the integrand is no longer positive but rather 
oscillating. 

Erdos' proof is done for the Poisson model. Recently, Klopp and Raikov 
|102j found a completely different approach based on the approximation by 
periodic potentials. Their proof works for alloy-type potentials. Moreover, 
the latter paper contains results for internal Lifshitz tails as well as for the 
case of an additional periodic background potential. 

We state Erdos' result: 

Theorem 3.11. If H u = Hq(B) + V w with a Poisson random potential 
and if the single site potential f > has compact support, then 



hxN(E) it 
lim K J = --. 3.54 

E\E InE B 

The reasoning we gave in the proof of Theorem l3.1Ul can be used to prove 
a lower bound for Theorem 13.111 as well. In fact, the bound suggests that 
the borderline between the two kinds of behavior is given by Gaussian single 
site potentials. This was actually proved in |t>4| and |39| . 

We have seen in this section that random Landau Hamiltonians may 
have unusual Lifshitz behavior compared to the nonmagnetic case. This 
originates in the fact that the lowest Landau band is flat, collapsing into 
one point. It is well known that this never happens for {B = and) periodic 
scalar potential (see |150j ). 

It is reasonable to expect that the flatness of the ground state Landau 
band is removed by an additional periodic scalar potential, at least generi- 
cally. If this is true, we would certainly expect "normal" Lifshitz behavior 
for such operators, including those with alloy-type potentials and q m \ n ^ 0. 
There are results in this direction in |102j . 
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We turn to the three-dimensional case with homogeneous magnetic field. 
We consider the Hamiltonian 

d 1 „ \ 2 ( . d 1 \ 2 d 2 



H (B) = i- -Bx 2 + i— + -B Xl 



dx\ 2 J \ dx2 2 J dxs 2 

In d = 3, the magnetic field itself introduces an anisotropy. The two space 
dimensions perpendicular to the magnetic field (37, £2) will be denoted by 
x±, the direction x% parallel to the field by xn. If we add an anisotropic 
Anderson (or Poisson) potential with a single site potential / obeying 

< /(*) < r: is b iz . ( 3 - 55 ) 



\X±\ ai + |X|||° 2 ~ ~ \x±\ ai + \X\ 

we will have Lifshitz behavior that resembles the results of Section 13.41 
except that in the _L-direction the behavior is always classical. 
Following the conventions of Section 13.41 we define 

2 1 

7i = — , 72 = — and 7 = 71 + 72- (3.56) 

If / has compact support in xj_-direction, we set 71 = 0. To have the 
potential well defined, we need 7 < 1. 

Theorem 3.12. The Lifshitz exponent r] for Hq{B) + V UJ in dimension 
d = 3 is given by 

+ max < - , >. (3.57) 



1-7 [21 

This theorem was proven for / with compact support in |178j . The case of 
72 > I is considered in |6()j and the mixed classical-quantum case is taken 
from |73] , 

We remark that the lower bound given above (in connection with Theo- 
rem EUJ can be used in this case as well. As a test function in this bound, 
we use 

*oOe±,Z||) =Vto(z±)0o(a;||) ( 3 - 58 ) 

where ipo is the (2-d) ground state (|3.50|) and (f>o(x\\) = L~ x l 2 is the ||-ground 
state. As before ^0 has to be cut down to zero near the boundary of the 
cube A^. This gives an error term of the order e~ L in the _L-direction and 
an error term of the order L~ 2 in the ||-direction. 

For further references about magnetic Lifshitz tails, we refer to fTT8] 
and to [TTfi] . 

There are also results on the integrated density of states and on Lif- 
shitz tails for random magnetic fields. We refer to the works of Ueki and 
Nakamura fTT2l ITTsl ITT41 IT381 IT39] and the references given there. 
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3.9. Percolation Models. Recently the integrated density of states 
for Laplacians on percolation graphs has been investigated. We consider 
bond percolation on Z rf with d > 2. This means that we remove bonds in 
the graph Z d independently with probability 1 — p for < p < 1. The 
resulting random graph is called the percolation graph and is denoted by Q p 
or simply by Q. 

For p small (p < p c ), all connected components ("clusters") of Q p are 
finite almost surely. For p > p c there is a unique infinite cluster (see, e.g., 

M) 

We consider the Laplacian on this graph which is a random operator, 
due to the randomness of the underlying graph. For a general discussion of 
percolation Hamiltonians, we refer to |176j . 

Actually, there are various Laplacians on Q due to different boundary 
conditions. To define these operators, we start with the adjacency operator 
Ag defined by the matrix elements with Ag(i,j) = 1 if \i — j\ = 1, i,j G Q 
and Ag = otherwise. For any i € Q, we let dg{i) denote the number of 
sites in Q to which i is connected. We denote by Dg the diagonal matrix 
with entries dg(i) on the diagonal. The Neumann Laplacian Lg is defined 
by 

Lg = Dg — Ag. (3.59) 
The Neumann Laplacian is the intrinsic Laplacian of the graph Q. In a 
sense it "ignores" the embedding of Q into Z d . The Dirichlet Laplacian Lg 
is defined by 

Lg = Dg + 2(2d- Dg) -Ag. (3.60) 
On its diagonal this operator counts the connections to neighboring sites 
in Q once and the (lost) connections to sites in Z d \ Q twice. It was an 
observation of Simon [T56] that the ab ove operator is a good analog of the 
Dirichlet Laplacian for subgraphs of Z d . See also |78j or |7()j for a discussion 
of boundary conditions for discrete Laplacians. 

Both for Lg and for Lg, a integrated density of states can be defined in 
analogy to Theorem 12.21 We call them Ng (E) and Ng(E), respectively. 

Since there are infinitely many clusters containing just one point, the 
Neumann Laplacian Lg has an eigenvalue of infinite multiplicity at the 
bottom E = of its spectrum. This causes the integrated density of states 
to jump at that energy. Hence 

^({0}) = <(0)-<(0-) = <(0) ^0. 

Theorem 3.13. (1) If p < p c , we have 



In -ln(A^(£)-A^(0)) 1 
lim — ^ — -L = --. (3.61) 



(2) If p > Pc, we have 



iim ln(^( E )-<(0)) =i 

E\o ln(E) 2 y ' 
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(3) For arbitrary < p < 1 , we have 

ln(-lnN°(E)) d 

lim V ~. = --. (3.63) 

E\o ln(E) 2 v ; 

Part (1) of this theorem, as well as part (3) for p < p c , was proven in 
|78j : part (2) and part (3) are taken from |133j . 

The perhaps surprising behavior of Ng (E) for p < p c is due to the fact 
that long one-dimensional chains dominate. They lead to small eigenvalues 
for the Neumann Laplacian since they "don't know" they are actually in 7L d . 
For the Dirichlet Laplacians these one-dimensional chains have rather high 
eigenvalues due to the additional diagonal term. In fact, they dominate the 
scene at E = 4d, the top of the spectrum. This can be seen by a symmetry 
argument. 

Part (2) of the theorem comes from the fact that for p > p c the leading 
behavior of Ng comes from the infinite cluster. In a sense, the infinite 
cluster looks rather cf-dimensional. The proof in [133] of this part relies on 
a celebrated paper by Barlow [5] . 

4. Regularity of the Integrated Density of States 

4.1. Introduction. In this final section we discuss regularity proper- 
ties of the integrated density of states. So far we have seen that the density 
of states measure v is a positive Borel measure with a distribution function 
N(E). 

Of course, the name integrated density of states suggests that v (resp. 
N) should have a density n(E) in the sense 

N(E) = [ n(E') dE', (4.1) 
Je'<e 

u([a,b}) = [ n(E') dE' . (4.2) 

J[a,b] 

We will, indeed, prove this for the Anderson model under certain assump- 
tions on Pq, the probability distribution of the random potential. The con- 
tinuous case is more complicated. For this case we will state only some key 
results and refer the reader to the literature. A good account is the re- 
view article by Veselic |177j . However, we remark that there are important 
developments after this survey; we mention especially the recent preprint 

m 

We will see that the integrated density of states does not always have 
a density. We will therefore also look at weaker regularity properties of 
N. We are especially interested in the question whether the function N is 
continuous, which is the same as i>({E}) = for all E. It turns out that 
this is always the case for the Anderson model. However, for the continuous 
case, no such result is known. 
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4.2. Continuity of the Integrated Density of States. The results 
of this section are valid for general ergodic operators on f(Z d ). Analo gous 
results for M. d have not been proven so far. The most general result in 
this context is the following theorem proved by Craig and Simon |25j in a 
somewhat stronger form (see Theorem 14 .3(1 . 

Theorem 4.1 ([221 ED])- Let {V u (n): n e Z d } be an ergodic stationary 
real valued random potential satisfying 

E[log(l + |VU0)|] <oo. (4.3) 

Then the integrated density of states of the Anderson operator h^ = ho + V^j 
is a continuous function, i.e., ^({A}) = for all A £ R. 

This result can be proven with an elementary argument by Delyon and 
Souillard |30L 12 7j using a kind of "unique continuation" property of discrete 
Schrodinger operators. We use this idea in the proof of the following lemma 
to show that no eigenspace can be sufficiently degenerated to produce a 
jump of the integrated density of states. 

Let us denote by /i#(-) the projection-valued spectral measure of H, i.e., 
Hh(A) = xa{H)'i in particular, ^#({A} ) is the projector onto the eigenspace 
of H with respect to A. 

Lemma 4.2. We have 

dim( X A L Ran^({A})) < C L^ 1 . 

Proof. The set 

A^ = < i 6 A^l max \i v \ = L or max = L — 1 

^ u=l,...,d v=l,...,d 

consists of the two outermost layers of Ax,. A solution u of h^u = Xu is 

(2) 

uniquely determined inside A^ by its values on A^ . So, the dimension of 
XA L (Ran/i/j^ ({A})) is & t most the number of points in . □ 
With the lemma, we can prove the theorem. 

PROOF of the Theorem. By Proposition 12.11 and Theorem 12.21 we 
have 

!/({A}) = lim {2L \ l)d tr {xa l Hh* ({A})) • (4.4) 

If fi is an orthonormal basis of xa l (Ran///^ ({A})) and gi an orthonormal 

basis of (xA i (Ran/i^({A})))" L we have, noting that x AL (Ran/z hij ({A})) is 
finite-dimensional, 



tr (xAi/^,({A})) = 22(fi,XA L fJ>h„(W)fi) 

<dim( X A L (Ran /U/l J{A}))) (4.5) 

<CL d ~ l (4.6) 
hence (|4.4f) converges to zero and ^({A}) = 0. □ 
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In |25| Craig and Simon prove a stronger result than just continuity: 

Theorem 4.3 (Craig-Simon). Under the assumptions of Theorem \4-l[ 
the integrated density of states is locally log-Holder continuous, e.g., for any 
positive r there exists a finite constant C r such that 

\N(X)-N(X')\<C r |lQg|A 1 _ A/|| (4-7) 

for |A| < r and \X - A'| < 1. 

The theorem of Craig and Simon is based on the Thouless formula for 
a strip in Z d . For the Thouless formula in the one-dimensional setting, we 
refer to the next subsection. 

4.3. Regularity of the DOS in Dimension One. There are very 
powerful techniques to study the Schrodinger equation and its discrete ana- 
log in dimension one which are unfortunately restricted to dimension one 
exclusively. Such techniques have been successfully applied to random op- 
erators as well to investigate the integrated density of states. 

At the heart of most of these techniques lies the reformulation of the 
eigenvalue equation of a second-order equation into an initial value problem 
for a system of first-order equations. Specifically, let us look at the eigenvalue 
equation 

hu = -u(n + 1) - u(n - 1) + (V u (n) - E)u(n) = 0. (4.8) 

We define 
and 

ME)= (V U n)-E -iy 

The function u(n) is a solution of ()4.8|) if and only if 

U(n + l) = A n+1 (E)U(n) 
for all n G Z. With the transfer matrix 

n 

$ n (E) : =n>(£) 

i=l 

the solution of (|4.8j) to the right initial condition 




can be expressed by 

U(n) = $ n (E)U(0). 

Similarly, it is possible to define the solution to the left. The spectral theory 
of the operator h^ is encoded in the matrices A n {E) (or Q n {E)). Note that 
these matrices belong to the group SX(2,]R). Thus it should not come as 
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a surprise that harmonic analysis on SL(2, R), in explicit or implicit form, 
plays a major role in the analysis of h w . 

The asymptotic behavior of the eigensolutions of h = ho + V is reflected 
by the Lyapunov exponent given by 

7 (£):= lim \n\\$ n (E)\\. (4.9) 

TV— >±oo 

This definition of the Lyapunov exponent is well defined by Fiirstenberg's 
theorem ( |46| . see also |27| ). The link between the Lyapunov exponent and 
the integrated density of states is expressed in the Thouless formula 

7 (E) = J \og\E - E'\ dN(E'). (4.10) 

This formula from physics ( [5611171] ^ was made rigorous in [3]. A simplified 
proof can be found in |26) or |27| . 

The Thouless formula and the obvious fact that "f(E) is non-negative 
imply that the integrated density of states TV" is log-Holder continuous (in 
one dimension). This was observed by Craig and Simon in |26j . The proof 
by the same authors of the multidimensional analog uses a version of the 
Thouless formula for strips in higher dimensions |25j . 

As mentioned in the previous section, these results hold for general er- 
godic potentials on Z d , they are (at least) close to optimal in this generality 
(see also our discussion below). 

However, if we assume that the Kj(n) are independent (and identically 
distributed) much more can be said about the regularity of the integrated 
density of states. In the next section, we discuss Lifshitz continuity in the 
multidimensional case under assumptions on the distribution Pq of the ran- 
dom variable K,(0). In the one-dimensional case we discuss here, we have a 
fairly complete picture about the regularity of N. For simplicity we assume 
that suppPo is compact. 

Definition 4.4. We call a function f onM Holder continuous of order 
a if 

\f(x)-f(y)\<C\x-y\ a . 

Theorem 4.5. Let h = ho + V be a discrete random Schrodinger oper- 
ator in dimension one. Assume V is a sequence of independent, identically 
distributed random variables, such that their distribution Pq has compact 
support. Then the integrated density of states N of h is Holder continuous 
of some order a > 0. 

This theorem is due to |114j . The simple case of the discrete Laplacian 
/io (i.e., V = 0) shows that one cannot expect more than Holder continuity 
in general. For ho the integrated density of states is merely Holder con- 
tinuous of order ^ at the band edges. The derivative of N, which is called 
the density of states, diverges at E = ±2. This behavior is known as van 
Hove singularities. The Lifshitz behavior of the integrated density of states 
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seems to suggest that N is smoother at the band edges for a truly random 
V. 

One can expect that more regularity for Pq implies more regularity for 
N. In fact, Simon and Taylor |161j proved that already a little regularity 
of Pq implies that N is C°°. To state their result, we need the following 
definition of a Sobolev space: 

Definition 4.6. We say that a function f G L p belongs to La if the 
Fourier transform f satisfies: There is a function g £ LP such that g(k) = 
(1 + \k\ 2 ) a / 2 f(k) is an TP-function. 

It is not hard to see that the characteristic function xi of a finite interval 
I is in L a for a < \ (see |161j ^. 

Theorem 4.7. Let h = ho + V be a discrete random Schrodinger oper- 
ator in dimension one. Assume V is a sequence of independent, identically 
distributed random variables. If the distribution Pq o/ 14,(0) has a density g 
with compact support and such that g G L a for some a > 0, then N E C°° . 

If we do not assume any regularity of Pq, Theorem 14.51 is the best result 
we can hope for. To see this, let us look at a Bernoulli distribution for Pq. 
We set 

P =p5 a + (l-p)5 b . (4.11) 

Halperin |55| proved, with some points of rigor clarified by Simon and Taylor 
QH], that: 

Theorem 4.8. Assume that the V w {n) are independent with identical 
distributions Pq as in 1)4. lip with < p < | . Then the integrated density of 
states is not Holder continuous of any order a larger than 

2|log(l-p)| 

CtQ = ; . 

arccosh(l + ^\b — a\) 

If N is Lifshitz continuous (Holder continuous of order a = 1), then iV 
has a bounded density n, i.e., N(E) = n(X) d\ and vice versa. If iV is 
Holder continuous of (strict) order a < 1, then N may still have a density 
which then has to be unbounded. So, Theorem 14.81 does not rule out that 
the density of states measure v is absolutely continuous. However, it is true 
that v has a singular continuous component if \b — a\ is large. 

Theorem 4.9. If Pq is Bernoulli (as in (|4.1ip ) and < p < 1, then the 
density of states measure v has a singular continuous component if \b — a\ 
is large. 

This theorem was proved in |14j following ideas from |161j . The paper 
|14j contains the first proof of Anderson localization for the one-dimensional 
Bernoulli model. The proof of a singular continuous component of iV is 
based on this knowledge. The article |122j continues these investigations 
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and proves that for \b — a\ large, the density of states measure is even purely 
singular continuous. 

For further results concerning the regularity of the DOS of random 
Schrodinger operators and its discrete analog in dimension one see \12\ \28\ 

mnimi9mi^i^i^i94iiT2miT7^ . 

As the operators with i.i.d. potentials, almost periodic operators belong 
to the class of ergodic operators. In a sense these two types of operators form 
the extreme cases within the class of ergodic operators. There are recent 
results on regularity of the integrated density of states for one-dimensional 
almost periodic operators by Goldstein and Schlag |51l 153) . We refer to the 
review |52j for details. 

4.4. The Wegner Estimates: Discrete Case. One motivation in 
physics to study the integrated density of states was the hope to use it as an 
indicator for different spectral types. The aim was to distinguish pure point 
spectrum and continuous spectrum at the mobility edge. In some sense the 
conjectures discussed were inconsistent. Some expected a divergent density 
of states at the mobility edge, while others assumed a vanishing density of 
states. 

In 1981 Wegner |179j put an end to this discussion by proving upper and 
lower bounds of the density of states in the discrete setting of the Anderson 
model. These estimates imply the density of states neither vanishes nor 
explodes at a mobility edge (or anywhere else). 

Wegner's result (more precisely, his upper bound) soon became a cor- 
ner piece in the proofs of Anderson localization by the multiscale analysis 
method and it still is. To formulate Wegner's result, we introduce boundary 
conditions on the lattice. The simplest boundary condition on the lattice is 
defining through its matrix elements: 

ht(i,j) = K(i,j) (4.12) 

whenever both i and j belong to A. For a discussion of boundary condition 
on £ 2 , see Simon |156j or the review |7()j . 
We define 

N A (E) = #{n\ E n (ht) <E} = tr (P^jj, (/£)) . (4.13) 

Theorem 4.10 (Wegner-estimate) . Suppose the measure P$ has a bounded 
density g, (i.e., Pq{A) = J A g(X)dX, \\g\\oc < oo), then 

E(N A (E + e)- N A (E - e)) < 2 IMU |A| e. (4.14) 

Remark 4.11. The assumption that Pq has a density cannot be dropped. 
We have seen at the end of Section \4-S\ that the integrated density of states 
has a singular continuous component if Pq is a (certain) Bernoulli distribu- 
tion. 

Before we discuss this estimate, we note two important consequences. The 
first concerns the regularity of the density of states. 
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Corollary 4.12. Under the assumption of Theorem \4-l(J[ the integrated 
density of states is absolutely continuous with a bounded density n(E). 

Proof. 

N(E + e) — N(E - e) = lim —M(N A (E + e) - N A (E - e)) 

|A|— >oo \A\ 

< Ce by Theorem EL"TU1 □ 

Thus N(E) = n(X) dX. We call n(X) the density of states. Some- 
times, N also is called the density of states which, we admit, is an abuse of 
language. The second consequence of Theorem I4.1UI is a key ingredient in 
proving Anderson localization. 

Corollary 4.13. Under the assumptions of Theorem \4.1C\ we have for 
any E and A, 

P(dist(.E, <T(h£)) <e)< Ce\A\. (4.15) 
Proof. By the Chebyshev inequality, we get 

P(dist(£, a{h^)) <e) = ¥(N A (E + e) - N A (E - e) > 1) 
<E(N A (E + e)-N A (E-e)) 
<Ce\A\ by Theorem I4~TU1 □ 

The first step in the proof of Theorem I4.1UI is to average the eigenvalues 
inside the interval (E — e, E + e) with respect to the random potential. To 
do this we consider the eigenvalues E n {h^) as functions of the arguments 
Vi = V u (i) with i E A, i.e., E n (h^) = E n (Vi,i & A). The resulting estimate 
is summarized in the following lemma. 

Lemma 4.14. With e > let g be a non- decreasing C°° -function with 
g(X) = 1 for A > e, g{\) = for X < -e and < g(X) < 1. Then 

r-E+2e f) 

N(ht E + e)- N(ht E-s)<J2 ^7 trg(^({V(}) - ^drj. 

Now we are in a position to interchange the expectation of the random 
potential and the energy integral. Since the random variables V w (i) are 
independent and have the common distribution dPo(Vi), the expectation E 
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is just integration with respect to the product of these distributions. Hence 



E(N A (E + e) - N A (E - e)) 

< E(]T f E+2£ A tr (e(h£M}) ~ v)dn)) 

r E+2e p a 

= / ^EII / ffW)%tr e (to)-,) 

JE-2e j £ M £ K J ^ OV i 

rE+2e r rb a 

<IMIoo/ d V J2 9{Vi)dVi dVj—trgihtm})- 

Je-2s f£ .J^. Jr Ja dVj 

rE+2e r 

^ ^e n a(Vi)dv t 

JE-2e je AieA,i& J 

{tr Q(ht({V^= b ) - V) ~ tr e(^({^ =a ) - V)} ■ 



Above, a, b are such that suppg C [a,b]. The notation {V^}^ means the 
family {V^}j g A with V% = Vi for i ^ j and V^- = b. The problem is now to 
estimate the trace difference. In the discrete context of the Anderson model, 
the variation of a potential value at one site is a rank one perturbation. 

Lemma 4.15. Let A be a self-adjoint operator bounded below with purely 
discrete spectrum Eq < E\ < . . . (where the eigenvalues are repeated accord- 
ing to multiplicity) . If B is a symmetric positive rank one operator, then 
A = A + B has eigenvalue E n with E n < E n < E n +\. 

Given the lemma, we now continue the proof of the theorem. We set 
A = H A ({Vi} v i =a ) and A = H A ({Vi} v J =b ). Obviously, their difference is a 
(positive) rank one operator 



tr g(A - n) - tr g(A - n) = ^2{g{E n - rj) - g(E n - r/)) 

< X)(e(^n+i - V) - Q{E n - r,)) 

< sup g(X) - g(fi) 

= 1. (4.16) 



We conclude the proof of Wegner's estimate by proving Lemma 14.151 
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PROOF. B = c\h >< h\ with c > 0, i.e., B <p = c{h,tp) h for some h. 
By the min-max principle (see |150j ). 

E n = sup inf (ip, Atp) + c \(<p, h)\ 2 

llvll=i 

< sup inf (99, j4</?} 

llvll=i 

< sup inf Ay?) 

IMI=i 

= £„.+!• (4.17) 

□ 

Wegner's estimate is intimately connected with a method called "spec- 
tral averaging." Roughly speaking, spectral averaging says that taking ex- 
pectation with respect to random parameters will make the spectral measure 
absolutely continuous. Here is a typical example which comes from the the- 
ory of rank one perturbations, see Simon |158j and references given there. 

To formulate this result, we take a bounded operator h on £ 2 (Z d ) and set 
h a = h+ a5j for any fixed j S Z d . Note that the multiplication operator Sj 
is a rank one perturbation of h. We denote by m a the (projection valued) 
spectral measure of h a and set (A a (A) = (Sj,m a (A) Sj). We obtain 

Theorem 4.16 (Spectral averaging). 

J dfi a (E) da = dE, (4.18) 



i.e., 



J (J f{E)d i x a {E) \ da = J f(E)dE 

for all integrable f . 

Wegner's estimate follows from the spectral averaging result (see |158p . 

Spectral averaging was introduced in the theory of random operators by 
Kotani |109) who used it to prove Anderson localization. He used random 
boundary conditions in dimension one, but soon Kotani's trick was also used 
to prove localization in higher-dimensional systems. However, Kotani was 
not the first to prove a spectral averaging formula. Such a formula was 
known in the Russian literature earlier, e.g., to Javrjan |65j who proved it 
for "random" boundary conditions. Spectral averaging also plays a promi- 
nent role for continuous Schrodinger operators; see |21j and references given 
there. 

4.5. Regularity in the Continuous Case. To prove a Wegner es- 
timate for the continuous case is considerably harder than for the discrete 
case. In fact, only the alloy-type model and a few other cases can be treated 
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so far (for these cases see [21] and |41| ). For the alloy- type model one can 
carry over Wegner's original proof ( |69j ). 

However, the finite rank estimate (Lemma I4.15JI cannot be transferred 
directly to the continuum (see |67j ). Thus a direct analog of Wegner's 
approach only gives 

E{N A (E + e) - N A (E - e)) < C|A| 2 e . (4.19) 

The estimate (|4.19|) obviously gives no information on the regularity of N. 
However, it suffices as input to multiscale analysis to prove Anderson local- 
ization (this was observed in |121| ). 

The first to prove a Wegner bound with the "correct" volume dependence 
were Kotani and Simon in jllOj . They required the single site potential / 
to be the characteristic function of the unit cube. Later Combes and Hislop 
|18j relaxed this condition. 

Meanwhile, there is a large number of results on "generalized" Wegner 
estimates of the form 

E(N A (E + e) - N A {E - e)) < C\A\ k e a . (4.20) 

For k = 1, they imply Holder continuity of the integrated density of states 
with Holder exponent a. 

We refer to the excellent survey by Veselic |177j on the subject which 
describes the development until 2004. In addition, we mention the papers 
HI EHl [22l E31 EH ED HH ITom [TOTI [T08] for further reading. 

Very recently, Combes, Hislop and Klopp |20j published a result which 
includes (and improves) virtually all previous results. 

Theorem 4.17. Let qi be independent random variables with a com- 
mon distribution Pq of compact support. Let f be a non-negative single site 
potential of compact support. Then 

(1) If Pq is Holder continuous with Holder exponent a, then N is 
Holder continuous with the same exponent. 

(2) If Pq is Lifshitz continuous, then N is Lifshitz continuous as well. 
In this case N has a bounded density. 

The authors of |20) actually prove a more general theorem allowing the 
random variables to be dependent and including a magnetic field. 

4.6. Beyond the Density of States: Level Statistics. One may 

look at the energy statistics of a disordered system on a smaller scale than 
we do for the integrated density of states. This subject is common in the 
theory of random matrices since the days of Wigner and Dyson (see, e.g., 
|124j or |29j ). but is still in its infancy for random Schrodinger operators. 

Suppose we have a random Schrodinger operator H which we restrict to a 
cube Al by appropriate boundary conditions. We call the resulting operator 
H Al . For a fixed cube A^, the operator H Al has roughly \ Al\ energy levels 
around an energy E. So we may say that the averaged level spacing near 
E is it I . We now look at the eigenvalues around E under a microscope 
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zooming the averaged level spacing to 1. The keyword is "unfolding of the 
spectrum." By this we mean we look at the measure 

b 



ti L ([a,b]) = #{n;E n (H AL ) G 



E + -—,E + 



(4.21) 



|Al|' 

It is reasonable to ask whether there is a limit of when L goes to 
infinity. Moreover, if such a limit exists, what are the properties of the limit 
measure /x? 

Molchanov |132j proved the existence of this limit for a one-dimensional 
model in the continuum. Minami |131j investigated the multidimensional 
discrete case in the regime of Anderson localization. Observe that we have 
localization for all energies for Molchanov's model. 

We discuss Minami's case (Molchanov's case being similar). In the fol- 
lowing we give merely a rough sketch of Minami's result, leaving out many 
details — even assumptions and precise statements. A complete discussion is 
far beyond the scope of this paper. We urge the reader to look at the paper 
|131j to get a complete picture. 

Let us suppose we have an integrated density of states ./V which has a 
bounded density n, i.e., N(E) = n(A) d\. For energies near the bottom 
of the spectrum, Minami shows that /^l is asymptotically a Poisson measure 
of intensity n(E). This implies that for a < b, 

V(a L ([a,b] = k)) « n( ^~ a) V "^- a >. (4.22) 

Moreover, the random variables hl(Ai), . . . ,/i^(^4 m ) are approximately in- 
dependent. 

In an informal (but provable) sense, this means that the eigenvalues 
near E look like independent random variables with a uniform distribution. 
Especially we have 

F(fi L ([a, b] = 1)) » n(E)(b - a) e -»(^)(&-«), (4.23) 

P(/^([a, b] = 2)) « n(E)2( 2 fe ~ a)2 e -^)( fe -) (4.24) 
and consequently one may prove 

Mi([-x/2,x/2])<2 » L ([-x/2,x/2])<l) « {n{E fjf )/2 ~n{E)x. 

This says in a rough way that the differences of energy levels near E have a 
probability density which is strictly positive near 0. In physics terminology, 
there is no level repulsion. 

One expects that this is not the case in the energy region of extended 
states. So far nobody has proven the existence of extended states for the 
Anderson model. A fortiori, nobody has proven level repulsion in this case. 
However, level repulsion is well established for the classical Wigner-Dyson 
ensembles of random matrix theory ( \124\ I29j ) . 
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